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The distribution of unicyclic components in a random graph is obtained analytically. The number 
of unicyclic components of a given size approaches a self-similar form in the vicinity of the gelation 
transition. At the gelation point, this distribution decays algebraically, Uk — (4fc)~^ for fc 3> 1. As 
a result, the total number of unicyclic components grows logarithmically with the system size. 
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Random graphs underly processes such as polymer- 
ization 0, percolation 0], and the formation of social 
networks P, f7| ■ Random graphs have been extensively 
studied, especially in theoretical computer science 0, Q • 
Special families of random graphs have been also exam- 
ined, e.g., planar random graphs appear in combinatorics 
0, l8| and in physics ^] . The basic framework for generic 
rando m g raphs naturally emerged in two different con- 
texts Flory [1 nil and Stockmayer [13 modeled 
a polymerization process in which monomers polymer- 
ize via binary chemical reactions until a giant polymer 
network, namely a gel, emerges. Erdos and Renyi stud- 
ied an equivalent process in which connected components 
emerge from ensembles of nodes that are linked sequen- 
tially and randomly in pairs |l3l |. 

Different methodologies have been employed to char- 
acterize random graphs. Kinetic theory, specifically, the 
rate equation approach, was used to obtain the size dis- 
tribution of components 0| . Using probability theory, a 
number of additional characteristics including in particu- 
lar the complexity of random graphs have been addressed 

In this study, we show that the rate equation approach 
is useful for studying the complexity of random graphs. 
Our main result asserts that Uk, the average number 
of unicyclic components of size fc in a random evolving 
graph, is given by 

C/.W^^f^e--^:^. (1) 

n— 

The unicyclic components size distribution becomes self- 
similar as the gelation transition is approached. At the 
gelation point {tg = 1), the size distribution develops an 
algebraic large-size tail: L/fc(l) ~ (4fc)^^. This implies 
that at the gelation point, the total number of unicyclic 
components depends logarithmically on the system size. 

The random graph evolves from N disconnected nodes 
as follows. At each step, two nodes are selected at ran- 
dom and a link is drawn between them. This linking 
process is repeated ad infinitum, leading to an ensemble 
of components, defined as maximally connected sets of 




FIG. 1: A random graph with TV = 10 nodes and 8 links. The 
two components on the left are trees, the two components on 
the right are unicyclic. 



nodes. We consider multi-graphs where the two selected 
nodes need not be different so that self-connections are 
allowed. 

Each component has a certain number of nodes and 
links. The number of nodes is the component size. The 
number of links minus the number of nodes is the compo- 
nent complexity [lj| . Up to a sign, the complexity is the 
Eulcr characteristic of the component — a topological in- 
variant. Trees have complexity —1, unicyclic components 
have complexity 0, bi-cyclic components have complexity 
1, etc (Fig.HJ. 

This linking process is treated dynamically. Initially, 
there are N components of size 1 and complexity — 1. 
Links are drawn between any two nodes with a uniform 
rate, set equal to 1/{2N) without loss of generality. It is 
useful to consider linking between different components 
and linking within the same component separately. Two 
different components of size and complexity (i, I) and 
(j, m) respectively merge with rate ij /(2N), symbolically 
represented by the reaction scheme 

{i,l) + ij,m)^ {i + j,l + m+l). (2) 

When linking involves two nodes in the same (fc, n) com- 
ponent, the process is 

(fc, n) -> (fc, n -f 1) (3) 
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and it occurs with rate fc^/(2iV). The average number of 
hnks at time t is Nt/2; the average number of sclf-hnks 
is smaller by a factor N, i.e. it is equal to t/2. 

Let the number of components of size k and complexity 
n be Nk,n- Following the dynamical rules ©-©j the 
quantities Nk,n change with time according to the bi- 
variate Smoluchowsky equation [TgI llTj 



obtained either using the Lagrange inversion formula 
or alternatively, by performing the direct calculation 
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The initial condition is Nk^n — N5kfl5n^-i. This rate 
equation implies that the number of trees, unicyclic com- 
ponents, and bi-cyclic components are proportional to A'^, 
7V°, and N^^, respectively. The random graph consists 
primarily of trees and unicyclic components, with more 
complex components being rare. Our goal is to derive 
the distribution of unicyclic components. Given the re- 
cursive nature of the rate equations, it requires the size 
distribution of trees, that to leading order in N equals 
the total size distribution. 

We first recapitulate the computation of the total size 
distribution. Let Ck — Nk/N be the size distribution of 



_(fc-l)! (fc-2)! 
The size distribution is therefore 0, IT^ 



(9) 



ah components (formally, Nk = J^i^kd)- Since linking tion - (27r) 



At time tg = I, the system undergoes a gelation tran- 
sition: a giant component that eventually engulfs the 
entire mass in the system emerges. Close to the gela- 
tion time, the size distribution attains the scaling be- 
havior Cfc(t) ~ kZ^^^^{k/k^,) with the typical size scale 
/c* — (1 — t)^^. This size scale diverges as the gela- 
tion point is approached. The underlying scaling func- 
-i/2^-5/2g-z/2 gxhibits an exponen- 



within a component affects its complexity but does not 
affect its size, the corresponding (first two) terms in the 
rate equation |0J are irrelevant and the total size distri- 
bution satisfies the nonlinear equation 
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tial large-size decay in the pre-gel regime while at the 
gelation time it develops an algebraic tail: Ck ~ fc~^/^. 
This power-law behavior allows us to estimate the size of 
the giant component, Ng ^ A^^/a using the extremal 
statistics criterion -^X]fe>Ar ^fe ~ 1 |2Clj| . The time when 
the giant component emerges in a finite system is esti- 



mated from Ng ^ {1 — tg) ^, i.e. 



1 -t„ - 



The initial condition is Cfc(O) = Sk,i- The size distribution 
is obtained via the generating function 



F{z,t) = J2kckit)e'' 



The average size distribution of unicyclic components 
Uk = (Nk^o) is coupled to the total size distribution c^. 
From (01) we find that Uk satisfies the linear inhomoge- 
neous equation 



k>l 



This generating function evolves according to 
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subject to the initial condition F(z,0) = e^. Let us re- 
write the derivatives through Jacobians: ^ ^ 9(f,z) 



and 
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Using the relation 
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recast the nonlinear equation for F{z, t) into a linear 
equation ^ = 1 — for z{F,t). Integration over time 
yields z{F,t) = (1 — F)t + InF. (The integration con- 
stant InF follows from the initial condition F{z, 0) = e^.) 
Exponentiating this equality yields an implicit relation 
satisfied by the generating function 



The initial condition is C/fc(0) = 0. The first term on the 
right-hand side of Hl()() plays the role of a source — it rep- 
resents the formation of unicyclic components from trees 
via the linking of two nodes within the same component. 
Such linking occurs with rate fc^/2iV (recall that multi- 
graphs where two nodes can be connected by more than 
one link are considered). The next two terms account 
for changes in the size of a unicyclic component due to 
mergers with different components. 

Consider the average number of unicyclic components 
U — Uk prior to the gelation time. Summing 
Eqs. (|10|l we find that U satisfies 
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dU 1 
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We can expand C, in terms of F, yet we are seeking the 
opposite: F = Y,kck{t)C,^ . The size distribution can be 



where M2{t) = J2k ^^^^(i) is the second moment of the 
size distribution. Using M2{t) = F'{z,t)\z=o and Eq. © 
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FIG. 2: The average number of unicyclic components at the 
gelation point versus the system size A'^. 



we get M2 = (1 — t) ^. The total number of unicyclic 
components is therefore 



U{t) = - In 
^ ' 2 1-t 



(12) 



The number of unicychc components diverges as the gela- 
tion point is approached. The total number of unicyclic 
components at the gelation point is obtained from the 
estimate 1 — ~ N~^^^; it diverges logarithmically with 
the system size N |2ll2a| 



(13) 



This number is much larger compared with the average 
number of self-connections that equals 1/2 at the gelation 
point. The approach to the asymptotic behavior H13|) is 
shown using Monte Carlo simulations (Fig. [5J . The data 
represent an average over 10^ independent realizations. 

To determine the entire distribution Uk, it is useful to 
manually solve for small k. The expressions 
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suggest the following structure kUk{t) = Akt^ e^'°*. Sub- 
stituting this form and the size distribution © into H1U|) 
we obtain a recurrence for the coefficients 



1 k'' 
Ak = 

2 fc! 
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To solve this recursion, we introduce the generating func- 
tions A{z) = Efe>i^fee'=^ and H{z) = 
thereby recasting 11411 into 



fc>l fc! 



A 



dH 

dz 



(15) 



The generating function H{z) = F{z -f 1, 1) is expressed 
via the known generating function F{z,t) and therefore, 
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Re- writing this equality ashiH — H = z and differentiat- 
ing with respect to z we obtain ^ = H/ (1 — H), which 
is then inserted into fTB|l to give A = ^H{1 — H)^^. 
Combining this equation together with Eq. (|16|l . the co- 
efficients are evaluated similar to the derivation of the 
size distribution c^: 
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Thus, the size distribution of the unicyclic components 
is obtained analytically. 
Given the scaling behavior of the size distribution, we 
investigate the unicyclic size distribution in the vicinity 
of the gelation point. To obtain the large-fc behavior, the 
size distribution is re-written as 
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When fc ^ 1, the product approaches an exponential, 
e-i(i+i)/2fe^ and the summation can be replaced by inte- 
gration 
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The leading large-fc asymptotics of the unicyclic com- 
ponents size distribution is therefore 



Uk{t) 




(18) 



In the vicinity of the gelation point, this distribution ap- 
proaches the scaling form Uk{t) ~ fc^^4'(fc/fc») with the 
same scaling variable as the one underlying the total size 
distribution, fc, = (1 — i)~^. The scaling function is also 
similar in form 



*(z) = (4z) 



(19) 



The unicyclic size distribution has exponential tails both 
above and below the gelation transition and it becomes 
algebraic at the gelation point 



Uk[tg) ^ (4fc)- 



when 



fc > 1. 



(20) 
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This result is of course consistent with the total number 
of unicyclic components: U{tg) ~ J2k<N^/^ Uk{tg) leads 
to Eq. 

We now probe the complexity of the giant compo- 
nent. In the post-gelation region the gel mass (mass 
not contained in finite components) g{t) is given by 
5=1- Efcfccfc(t) - 1 - F{z,t%^a- Using Eq. ©, 
the gel mass obeys g — 1 — e^^* for t > 1. The com- 
plexity of the giant component increases due to linking 
processes involving its internal nodes. The total number 
of nodes in the giant component is Ng. The linking rate 
is therefore (2iV)~^ x (Ng)"^, so the complexity of the 
giant component is 

C{t)^- dt'g^it'). (21) 

Just above the gelation transition, the gel mass increases 
linearly, g ~2{t-l), and therefore, C{t) ~ | N{t - if 
as t [ 1. In the large time limit there are a few 
monomers apart from the giant component. The total 
number of links is Nt/2 and the total number of nodes 
is N — Ni = N{1 — ci), so the complexity of the giant 
component grows according to N{t/2 — 1-1- e~*) when 
t — > oo. This result also follows from the general formula 
(|21|l and the asymptotic behavior of the gel mass g{t). 

In this study, we considered size and complexity char- 
acteristics obtained via an average over infinitely many 



realizations of the linking process. For some quanti- 
ties, averaging is irrelevant in the thermodynamic limit 
— > oo. For instance, the tree distribution is an exten- 
sive random quantity {(Nk -i) = Nck) with fluctuations 
of the order of iV^/^ |2l,l3|, so relative fluctuations de- 
crease with the systems size according to N^^^^. Thus, 
for trees the average distribution Nck well represents the 
outcome of a single realization of the random evolving 
graph. For unicyclic components, fluctuations are of the 
same order as the average and an analytical computation 
of the correlation function {Ni,o{t)Nj^o{t)) is a challeng- 
ing open problem. This correlation function is required 
for the determination of the average distribution of bi- 
cyclic components, so the naive form of the rate equa- 
tion Q is inadequate for describing bi-cyclic (and more 
complex) components. 

In conclusion, we have obtained the size distribution 
of unicyclic components in a random graph. Overall, the 
unicyclic size distribution has similar properties as the 
total size distribution; It becomes self-similar near the 
gelation point, and generally, it has an exponential tail. 
Precisely at the gelation point the distribution has an 
algebraic tail. The main difference with the size distri- 
bution is in the value of the power-law exponent itself. 
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